Interfacial effects on the polarization of BiFeO^ films 
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Abstract 

By considering an interfacial layer between the electrode and the BiFeOz{BFO) layer, the 
polarization and the hysteresis behavior of BFO film are simulated. It is found that the non- 
ferroelectric interface will increase the coercive field, and remarkably suppress the polarization of 
the ultrathin film under low applied fields. Due to the competition between the interfacial effect 
and the internal compressive stress, the maximum polarization on the P-E loop of a BFO film can 
be independent on the film thickness under an adequate applied field. 
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I. INTRODUCTION 



Recently, the perovskite-type oxides, which display ferroelectric and magnetic properties, 
have attracted considerable interest for their potential applications in multifunctional mag- 
netoelectric devices! 111. It has been shown that BFO is a good candidate which exhibits 

a 

the coexistence of ferroelectric and antiferromagnetic orders with high Curie temperature [2| 
and Neel temperature {3]. Numerous studies have been performed on BFO samples and 
especially, more recently, on thin films. Generally, bottom and top electrodes are fabricated 
to form a capacitor configuration for ferroelectric measurements. The physical origin of the 
interfacial layer has been proposed to be process-induced domain evolution or the effect of 
lattice misfit between the electrode and surface of the film. The interfaces in BFO film 
capacitors have already been verified by using X-ray diffraction and scanning transmission 
electron _pe high-angle annn.ar darfc-field imagingQ, and the effect of the interface 
on the dielectric behavior of BFO film capacitor was also investigated in experiment |5|]. 
For epitaxial system, the substrate-induced strains have been proven to have great influ- 
ence on the ferroelectric properties of ferroelectric material 6j. The previous studies show 
that substrate-induced strains can lead to a substantial increase in spontaneous polariza- 
tion for BFO films The mechanical substrate effects and thickness dependence 
of ferroelectric properties in epitaxial BFO films are investigated within the framework of 



Landau-Devonshire theory. Although 
has been reported in experiment 



1 the interface between the electrode and the BFO layer 



101 ]. there is no theoretical investigation of the interfa- 



cial effect on the polarization of BFO film. In this paper, the interfaces between electrode 
and BFO layer are considered with different film thickness. It is found that the interface 
will have great influence on the polarization and the hysteresis behavior of ultrathin BFO 
films. Due to the interfacial effect, the maximum polarization on the P-E loop of thin film 
can be lower(or higher) than that of thickness film, and the maximum polarization can also 
be thickness-independent through the competition between the interfacial effect and the 
compressive stress under an adequate applied field. 
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II. MODEL AND THEORY 



Here, we consider a BFO film grown on SrTiOs substrate with a bottom SrRuO^ elec- 



trode and a top Pt electrode 
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121 ] . It is known that there is a clear and sharp interface 



between BFO layer and the bottom SrRuO^ electrode, and a rough, inter-diffused interface 
between BFO layer and the top Pt electrode, which can be confirmed by the cross-section 



HR-TEM image 



10] . It is reasonable that the interface between BFO and the bottom 



electrode can be neglected for simplicity. Therefore, we simulate the BFO film by a ferro- 
electric layer with a top non-switching interface, and investigate the effects of the interface 
on the P-E loop with different film thickness. The electric displacement across BFO layer 
is given by Df = ejEf + Pf, where e/, Ef and Pf denote the permittivity, the local electric 
field and the spontaneous polarization of BFO layer, respectively. The electric displacement 
across the interface is D in = ei n E in , where ei n and E in stand for the permittivity and the 
local electric field of interface. The applied field E ex can be written as 



E P 



)Ef + v in E in , 



(1) 



where Vi n is the thickness ratio of interface. The continuity of the total current density 
across the film can be expressed as 



OP, 



J = cr in E in + e in E in = ctjEj + (e/ + —J-)Ef 



dE f 



(2) 



Combining Eq. (1) with Eq. (2), we obtain 



<r in E ex + e in E ex = [(1 - v)a in + va f ]E f 
+ [v(e f + X f) + e m (l-v)}E f , 



(3) 



where a in and aj are the conductivities of top interface and BFO layer, respectively, %f is 
the electric susceptibility of BFO layer, which equals dPf/dEf 13] 



Xf = P f {25cosh 2 [(E f -E c )/25]}-\ 
5 = E c [log(- +Pfr/Pfs ^ 



(4) 



1 - P fr / P fs 

where Pf r , Pf s and E c denote the remanent polarization, saturate polarization and coer- 
cive field of BFO layer, respectively, which can be determined through Landau-Devonshire 
theory [ 



The BFO films with different thickness are supposed to grow on (OOl)-oriented SrTiO^ 
substrate. A tetragonal structure of BFO thin film based on cubic perovskite structure with 
symmetry lowered to P4mm(neglect the small monoclinic distortion of about 0.5° of c axis) is 



suggested with an elongated c axisflj]. The polarization for such a single phase is along with 
c axis. The thermodynamic potential of pseudo cubic BFO film is expressed as the function 
of magnetization M/, polarization Pf, temperature T and misfit strain u m = (a s — ao)/a s , 
where a s is the in-plane lattice parameter for substrate and ao is the equivalent cubic lattice 
constant of free standing film. We do not include the magnetic effect and the coupling 
between magnetic order and ferroelectric order, owing to a relative weak magnetic order 
compared to the stronger ferroelectric one. The thermodynamic potential of BFO layer can 
be written as following, 

u 2 

F(P f , T) = alP] + a* n P} + a lu Pf - E f P f + - f- , (5) 

20 E (O e ) 2 p 

where a* = «i + u m s , ot*\ = a u + gr+fe ' Q% * s ^ e e l ec t ros tricitive coefficient, Sij 
is the elastic compliances, and the dielectric stiffness coefficient ai — (T — Tc)/2eoC. The 
equilibrium state polarization Pf under a local field Ef, Pf r and Pf s in Eq. (4) can be 
numerically resolved from Eq. (5) by dF/dPf = 0, 

2a{P f + 4a* n Pf + 6a lu p 5 f -E f = 0. (6) 

In order to introduce the misfit strain into the epitaxial BFO film, an effective substrate 
lattice parameter a s is defined as a s = p ^ s ^ +l , which is used to calculate the misfit strain u m 
instead of actual substrate lattice parameter a s , and u m = as ^~^)^ ■ P * s the equilibrium 
linear misfit dislocation density at the deposition temperature T G . p = u ™^} (1 — ), where 
h p is the critical film thickness corresponding to the generation of misfit dislocation. All the 



coefficients in Landau expansion used in this letter are given in Ref. 15]. 



III. NUMERICAL RESULTS AND DISCUSSIONS 

From Eq. (6), we calculate the spontaneous polarization Pf without considering the in- 
terface. The hysteresis behavior of the single BFO layer is shown in Fig. 1, where solid line 
stands for the hysteresis loop of the 77-nm-thick film, and dotted line for the 960-nm-thick 
film. The remanent polarization of the 960nm BFO film is 64/iC/cm 2 , which is close to the 
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polarization of single crystal[l^]. We can see that the spontaneous polarization increases as 
the thickness decreases, and our results are qualitatively consistent with previous theoretical 
works p], |8[. It is obviously that the compressive stress between BFO layer and substrate 
SrTiOz will remarkable increase the polarization, and the spontaneous polarization Pf will 
increase with the decrease of thickness. However, as observed by AFM[4] and simulated 
by the capacitor configuration [5], there really exists an interface between top electrode and 
BFO layer. So, the measured polarization is not the polarization of a single ferroelectric 
layer. In the following, an interface between the top Pt electrode and BFO layer is in- 
troduced, and two BFO film samples of thickness 77nm and 960nm are selected in our 
simulation [3]. The dielectric constant and the conductivity of the BFO layer are selected 
as 100 1^1 and 1 x 10~ 10 fi -1 cm -1 11], respectively. Due to the interface is formed between 
the metal material and BFO layer, we assume its conductivity and dielectric constant take 
the values of 1 x 10 _9 fi _1 cm _1 [5[ and 20. The thickness of interface is introduced by taking 
into account the thickness ratio of interface to BFO layer. The ratio v will increase as the 
decrease of film thickness, and the ratios v are taken as 0.02 and 0.002 for the two samples. 

For a given sinusoidal field E ex , the effective field Ef can be calculated from Eq. (3) as a 
function of time t. Then the polarization Pf at the mth increment in time is calculated by 
using Pf m = Pf m _ 1 + {Ef m — Ef m _^)j¥L\ m _ x . The simulated hysteresis loops corresponding 
to different applied fields are shown in Fig. 2, the frequency of the applied field in our 
simulation is 2kHz. The solid lines in Fig. 2 denote the loops of 77-nm-thick film, and the 
dotted lines stand for the loops of 960-nm-thick film. From Fig. 2(a) to (d), the maximum 
field are selected as 1.0 MV/cm, 1.25 MV/cm, 1.6 MV/cm and 2.5 MV/cm, respectively. By 
taking into account the interface, the local field Ef will decrease as the decrease of thickness 
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under the same applied field. For the 960-nm-thick film, due to the small proportion of 
the interface component, the local field Ef is close to the applied field, and the P-E loop is 
already saturated in Fig. 2(a), whereas the P-E loop for 77-nm-thick film does not saturated, 
and the maximum polarization on the P-E loop is smaller than the saturated polarization 
of 960-nm-thick film. One can see, from Fig. 1 and Fig. 2(a), that the interface will great 
suppress the polarization of the ultrathin film, and the interfacial effect of an ultrathin film 
is more important than the internal stress effect under a low applied field. For a selected 
applied field E max =1.25 MV/cm in Fig. 2(b), the maximum polarization on the P-E loop 
is independent on the thickness of the film, which is the consequence of the competition 
between the interfacial effect and the internal compressive stress effect. In Fig. 2(c), the 
maximum polarization of the ultrathin film is larger than that of the thick film, which 
had also been observed in experiment (See Fig. 4(d) in Ref. 18|), and it shows that the 
interfacial effect is less dominant for the ultrathin film under a high applied field. As further 
increasing the maximum field, the polarizations for the two samples are both saturated, and 
the saturated polarization of the 77-nm-thick film is larger than that of the 960-nm-thick 
film. In addition, the P-E loops for materials such as BFO film are often dominated by 
leakage current, and the saturated P-E curve in Fig. 2(d) for an ultrathin film can only be 
observed for samples with high resistivity and high breakdown field. 

In summary, the epitaxial Pt/BFO/SrRuO^/ SrTiO^, film is simulated by a multilayer 
model which takes the effect of interface. It is found that the existence of interfacial layer 
will increase the coercive field, and change the shape of the P-E loop, especially for an 
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ultrathin film. The P-E loop will become more slant for the ultrathin film, which is resulted 
from the decreasing of the effective field in BFO layer. For different maximum fields, the 
maximum polarization of the ultrathin film can be smaller(or larger) than that of the thick 
film. For an adequate applied field, the maximum polarization on the P-E loop could be 
thickness-independent, and the origin is caused by the competition between the interfacial 
effect and the internal compressive stress. It is concluded that the effect of interface between 
the electrode and BFO layer has great influence on the polarization of BFO film. When 
measuring the saturated P-E loops for multiferroic ultrathin films with metal contacts, one 
should apply a field large enough to saturate the polarization. 
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Figure Captions 

FIG1: Polarization of BFO layer as a function of the applied field without considering 
the interface. 

FIG2: P-E loop of BFO film with interface for different applied fields. Solid lines corre- 
spond to the 77-nm-thick film, and the dotted lines correspond to the 960-nm-thick film. 
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